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Abstract 

We examine certain issues related to the universality of the SU(2) lattice 
gauge theory at non-zero temperatures. Using Monte Carlo simulations and 
strong coupling expansions, we study the behavior of the deconfinement tran- 
sition in an extended coupling plane (/?, around the tricritical point where 
the deconfinement transition changes from second to first order. Our numer- 
ical results on A'^^ =2,4,6,8 lattices show that the tricritical point first moves 
down towards the Wilson axis and and then moves slowly upwards, if at all, 
as the lattice spacing is reduced. Lattices with very large A''^ seem to be 
therefore necessary for the mixed action to exhibit the critical exponents of 



*Email:gavai@mayur.tifr.res.in 

^Permanent address: S. N. Bose National Centre for Basic Sciences, JD Block, Sector III, Slat 
Lake, Calcutta 700064, India. Email:manu@hpth. difi.unipi.it 



1 



the three dimensional Ising model for positive values of the adjoint coupling. 
PACS numbers: 11.15Ha 
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1. INTRODUCTION 



Confinement of the non-abelian color degrees of freedom has been a challenging problem 
ever since gauge theories were formulated for quark-gluon interactions. In the beginning, our 
understanding of such gauge theories followed mainly from perturbation theory. Due to the 
property of asymptotic freedom, the perturbative approximations are valid for short distance 
phenomenon but they are inadequate to explore long distance physics like confinement of 
quarks. A strong evidence for, and a much better insight of, the color confinement mechanism 
in non-abelian gauge theories has been provided by analytical computations and Monte 
Carlo simulations of quantum field theories with a non-perturbative lattice regularisation. 
The simplest of such theory is described by the SU{2) Wilson action IQ. It was expected 
that non-abelian gauge theories in general do not have any phase transitions separating 
strong and weak coupling regimes. Therefore confinement, explicitly shown on the lattice in 
the strong coupling region, should persist also in the continuum limit. However, it was later 
found that certain non-abelian lattice gauge theories (e.g. 5*0(3), S'f/(4), SU{5)), with 
Wilson form of action undergo bulk transitions separating strong confining region from the 
weak coupling region where the continuum limit of the theory exists. Bhanot and Creutz 
P], extending the form of the action proposed by Wilson, showed that this apparent loss 
of confinement can be attributed to lattice artifacts associated with the choice of action, 
namely the so-called bulk phase transitions. 

Subsequent to the work of Bhanot and Creutz to characterize the bulk transitions in the 
extended coupling plane, important reasons for further exploration of this action in the past 
have been to study the basic mechanism of confinement itself ^ and to find out the origin 
of these bulk transitions 0-0. It has been a common folklore that the abrupt change from 
the strong coupling region to the scaling region for the Wilson action is due to the proximity 
of the critical point at the end of the first order line CD in Fig. 1, where the phase diagram 
obtained in Ref. |^ is shown, and that a "smoother" continuum limit may be obtained 
by going to negative (3a, the additional coupling for this action. Due to the theoretical 
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expectations of the role of S0{3){SU{N)/Zn) monopoles in SU{2){SU{N)) confinement [§, 
the above model is tailor made to study the interplay of these topological degrees of freedom 
and their role in confinement between SU{2){SU{N)) and S0{3){SU{N)/Zm) lattice gauge 
theories. In fact, the plaquette susceptibihty peak in the cross-over region in SU{2) lattice 
gauge theory and the bulk transition line BCD have been attributed to the underlying 
5*0(3) theory and its first order transition P,p!o|. These issues can be also analyzed and 



tested by exploring the mixed action at non-zero temperatures and will be further discussed 
after presentation of our results. 

The rich phase diagram associated with the mixed action, shown in Fig. 1 by solid 
and dashed lines, was established mainly by Monte Carlo simulations on relatively small 
lattices ||^ (4^-5"^) with periodic boundary conditions. Since these small lattices were also 
at finite temperature, the phase diagram is incomplete in the absence of the deconfinement 
transition line. Along the /?a = 0.0 axis, several finite temperature investigations have shown 
the presence of a second order deconfinement phase transition. Its critical temperature has 
been shown [|rT| to exhibit asymptotic scaling and its critical exponents have been shown 



to be in very good agreement with those of the three dimensional Ising model. Effective 
field theory arguments for the order parameter were used by Svetitsky and Yaffe [|r^ to 
conjecture the finite temperature SU(2) gauge theory and the three dimensional Ising model 
to be in the same universality class. The verification of this universality conjecture thus 
strengthened our analytical understanding of the deconfinement phase transition. Our work 



l^ , p!5[| on the extended action at non-zero temperatures began with the motivation to locate 
the line of deconfinement transition in the couplings plane Our simulations yielded 

the following surprising results: 

a] The transition remained second order in agreement with the universality conjectured 
exponents up to /3a ~ 1.0 but it became definitely first order for large enough /5a (> 
1.4). 

b] There was no evidence of a second separate transition at larger Pa, as would be sug- 



gested by the claim of Ref. 0] of a bulk transition there. For A^,- = 4 lattices the line 
of deconfinement transition was coincident with the line of bulk transitions of Ref. 
but for Nr = 2 there were no symptoms of any transition at those locations. The line 
of deconfinement phase transition, on the other hand, did move to smaller (3 for all 
13a, as Nt changed from 4 to 2. 

While the details of our analysis and results can be found in the works cited above, the 
key findings which lead us to these conclusions were following: 

a] The deconfinement order parameter, (L) (see next section for definition), acquired 
nonzero large value at the only transition found on all lattices studied (i.e A^t-=2,4,6,8) 
and showed clear co-existence of both phases at the transition point for larger 

b] The same critical exponent which established the transition to be in the Ising model 
universality class for 0.0 < [3a < 1-0 became equal to the space dimensionality, as a 
first order deconfinement phase transition would have, for larger 13a- 

c] The plaquette susceptibility showed a decrease at = 1.1 when the lattice 4- volume 
was increased by a factor of 16; it should diverge, i.e., increase 16-fold, if there were a 
first order bulk phase transition at (3a = l-l- 



The plaquette susceptibility results above are very similar to those of Ref. who too 
found a decrease in it while increasing the lattice volume by a factor of 16. On increasing 
the lattice size further, no further finite size dependence was found, leading to a conclusion 
that the finite size effects on smaller lattices are due to finite temperature effects. While 
larger lattices will be needed in our case too to see if a similar conclusion is reached, it has 
to be emphasized that conclusions based on finite size scaling usually do assume that the 
lattices are large enough for the scaling to set in. Thus a distinguishing feature between 
the bulk and deconfinement phase transitions, i. e. the finite size scaling behavior of the 
coupling at which the transition takes place with the temporal size of the lattice which 
leads one to expect the bulk transitions to move much less compared to the deconfinement 



transitions, is not necessarily useful here since it is not clear how big lattices are needed for 
this behavior to set in at various values of Pa- We have therefore relied heavily on the order 
parameter (L) to label a transition as a deconfinement phase transition, as mentioned in a] 
above. 

Recently, the above surprising results showing the change in the order of the deconfine- 



ment transition and the absence of the bulk transition were confirmed ||16| for another variant 
of the SU(2) action with a Villain form for the adjoint S0(3) part. We will later comment 
more on the above action. Taken together, these results pose many questions about the con- 
tinuum limit of the deconfinement phase transition and about the existence of separate bulk 
phase transitions. The foremost amongst them is about an apparent qualitative violation of 
the universality^, since an apparently irrelevant coupling seems to change the order of the 
deconfinement phase transition. The early simulations of the SU{2) lattice gauge theory 
are known to have yielded quantitative violations of universality [|18|. However, attributing 
them to the ignored higher orders in g"^, it has been shown [p!9|-|2T[] that dimensionless ratios 
of physical quantities have much weaker such violations. Indeed, one can hope that these 
violations will smoothly disappear under the error bars of the simulations as the cut-off 
becomes smaller. This will obviously not be the case for any qualitative violations of uni- 
versality. Of course, the region of couplings where universal results are obtained may have 
such an irregular shape that still larger lattices are needed to obtain universal results. It is 
not clear in that case, however, what the universal result would be. Clearly, if a universal 
result exists in the a ^ limit, then the tricritical point T, where the deconfinement phase 
transition changes order, must not appear on any renormalisation group flow lines to the 
Qu = 0.0 critical point and must thus be invisible in the continuum limit. If the point T 
moves up to large positive (3 a with increasing A^^ (and decreasing lattice spacing a), then the 



^Note that this universality, which results from the freedom of choice of the lattice type and 
action, is different from the finite temperature universality of critical exponents discussed earlier. 
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SU(2) deconfinement phase transition could still be of second order in continuum limit with 
Ising model exponents. This would be so irrespective of Pa used for simulations. If, on the 
other hand, the point T moves to large negative Pa, the universality with Ising model will 
be lost and the transition will be first order, again irrespective of Pa used in simulations. 

It may be argued that the presence of a line of bulk phase transitions and its end point 
will strongly modify the approach to continuum limit and thus large lattices are mandatory 
for seeing the universal physics at large Pa- It needs to be noted therefore that the bulk 
line in question was established only in numerical simulations on small, 4^-5^ lattices. A 



recent simulation p2| at Pa = 1-25 on larger A^^ lattices, with = 6, 8, 10, 12 and 16, 
found 1) a linear decrease in the average discontinuity in plaquette, AP, with A^ and 2) a 
plaquette susceptibility exponent of 2.09± 0.31 in contrast to the expected value of 4 for 
a first order bulk phase transition. This suggests that the end point of the bulk line is at 
Pa > 1.25. This is explicitly shown in Fig. 1 by drawing solid and dashed bulk lines above 



and below Pa = 1.25 respectively. While the result of Ref. ||2^ does explain the above 
mentioned results on plaquette susceptibility, the mystery of the apparent coincidence of 
the two different transition lines still remains for larger Pa- As we will show below, the 
deconfinement phase transition for A^- = 4 lattices, as identified by the order parameter, 
turns first order already at Pa = 1.25, suggesting that the bulk line or its end point 
are unlikely sources of this change. 

In this paper, we address the issue of the trajectory of point T with decreasing lattice 
spacing a, after defining in the next section the action we investigate and the observables 
we use along with their scaling laws. A simple strong coupling calculation is presented in 
Section 3, which suggests that the point T moves up in the plane to infinity. However, in 
our numerical simulations, described in Section 4, we find that it moves down on going from 
At- =2 to 4. On increasing A^,- further to 6 and then 8, we observe a very small upward 
movement by comparing the relative shapes of the Polyakov loop histograms. The last 
section contains a brief summary of our results and their discussion. 
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2. THE MODEL AND THE OBSERVABLES 



The lattice action is constrained only by a) the gauge invariance and b) the limit of zero 
lattice spacing which must coincide with the continuum form of the action. Infinitely many 
different forms satisfying these criteria can be written down. Bhanot and Creutz extended 
the Wilson action to a form described by the action, 

^ = E (/3 (l - IttfUp^ +Pa(i- I^taUp^^ . (1) 

Here Up denotes the directed product of the basic link variables which describe the gauge 
fields, f/^(x), around an elementary plaquette P. F and A denote that the respective traces 
are evaluated in fundamental and adjoint representations respectively. Comparing the naive 
classical continuum limit of eq. (|1|) with the standard SU{2) Yang-Mills action, one obtains 

1 = ^ + ^. (2) 

Here Qu is the bare coupling constant of the continuum theory. Introducing another coupling 
9, defined by tan6' = Pa/P, the asymptotic scaling relation [^] for this action is 

1 



- A(^) 

where 



exp 



Pogl]^ , (3) 



A(0) _ 57r2 6tan^ 
°^A(^ ~ 11 (3 + 8tan^) ' 

Here /3o and f3i are the usual first two universal coefficients of the /3 function for the SU{2) 
gauge theory: they do not depend on 9. 

One sees clearly from the equations above that the introduction of a non-zero (3 a, leads 
merely to a different and a correspondingly different value for the scale A{9). However, 
each of these theories, including the usual Wilson theory for /Sa = 0.0 flow to the same 
critical flxed point, = 0, in the continuum limit and has the same scaling behavior near 
the critical point. The different forms of action, obtained by varying Pa, are simply related 

8 



by a redefinition of coupling constant and tlie intrinsic scale A and yield the same universal 
continuum physics. Numerical investigations for different (3^ thus constitute a necessary 
check of the finite cut-off effects in the non-perturbative results obtained for Pa = 0.0, i.e., 
the Wilson action. 

Bhanot and Creutz [Q] found that the lattice theory defined by the extended action has a 
rich phase structure (Fig. 1). Along the (3 = axis it describes the 5*0(3) model which has a 
first order phase transition at PJ^^ ~ 2.5. At /?a = oo it describes the Z2 lattice gauge theory 
again with a first order phase transition at P"^"^ = |£n(l + \/2) ^ 0.44 p4|. Ref. found 



that these first order transitions extend into the {P,Pa) plane, ending at an apparent critical 
point located at (1.5,0.9). These transition lines are shown in Fig. 1 by continuous lines. 



Using finite size scaling, Ref. p2[ has recently shown that the critical endpoint must have 
Pa > 1.25. More simulations on larger lattices will be required to determine the endpoint 
precisely. The qualitative aspects of this phase diagram were also reproduced by mean field 
theory and large N and strong coupling ^ expansions. 

Simulations of the mixed action above at finite temperature are made on asymmetric A*"^ x 
Nr lattices, with periodic boundary conditions in the (shorter) r-direction. The partition 
function at finite temperature is given by, 

Z = [ YldU^ix) exp{-S) . (5) 

The order parameter for the deconfinement transition is the Polyakov loop pSf defined by 



1 Nr 

L{n) = -Trl[Uo{n,T), (6) 

^ T=l 

Here f/o(n, r) is the time-like link at the lattice site {n,T). Due to periodic boundary 
condition in the time-like direction at finite temperature the action of eq. (|T]) has a Z2 
invariance corresponding to the center of the gauge group. Defining this symmetry to be 

Uo{n,To) — * zUo{n,To) Vn, tq : fixed , and z E Z2 , (7) 

one sees that under its transformation the Polyakov loop changes by 



zL , (8) 

while the action in eq. (|l]) remains unchanged. 

A non- vanishing value for (L), with respect to the partition function in eq. signals a 
spontaneous break-down of the global Z2 symmetry. (L) is also an order parameter for the 
deconfinement phase transition, as it (or equivalently its average value L = -^^^ L{n)) can 

n 

also be shown to be a measure of the free energy of an isolated free quark In order to 
monitor the critical behavior of the deconfinement transition, we also define the Polyakov 
loop susceptibility: 

XN. = Nl{{L') - {Lf) . (9) 

In the thermodynamic limit, a second order transition is characterized by the following 
critical exponents: 

(L) oc |T - for T ^ (10) 

Xoc|T-Te|-^ for T^Te (11) 
^oc|T-T,|^'' for T^T, . (12) 

Here ^ is the correlation length corresponding to the Polyakov loop correlations and 
13 ~ 0.325, 7 ~ 1.24 and v 0.63 are the Ising model exponents, assuming the universality 
conjecture to be true. The best determination of these exponents for the SU{2) lattice gauge 
theory was made by using the finite size scaling theory according to which, the 



peak of the L-susceptibility on a lattice of spatial extent N„ is expected to grow like 

xnT^n: , (13) 

where = 7/z/ = 1.97. If the phase transition were to be of first order instead, then 



one expects the exponent = 3, corresponding to the dimensionality of the space In 
addition, of course, the order parameter is expected to exhibit a sharp, or even discontinuous, 
jump and the corresponding probability distribution should show a double (multi) peak 
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structure. For (3a = 0, the universality prediction was verified by Monte Carlo simulation 
by Engels et. al. |T^, who found u = 1.93 ± 0.03, whereas we found u = 3.25 ± 0.24 for 
/3a = 1.4 on = 2 lattices. 



3. STRONG COUPLING 



Before turning to the results of our simulations to determine ui and to locate the tricritical 
point, it may be an instructive exercise to find out what hints the strong coupling expansion 
can provide. Such expansions for the free energy [^,^ and string tension have been 



used in the past to study SU{N) deconfinement transition for the Wilson action. The 
basic strategy is to obtain an effective potential for the order parameter L, by expanding 
the partition function in powers of the inverse coupling constant (s) and integrating out the 
spatial links . Due to the Z2 symmetry of the theory, the Landau- Ginzburg effective action 
is an even polynomial in the Polyakov loop for the SU{2) theory. To lowest order: 



Seff = E log (1 - in)) - 4 ^j:L{n)L{n + 



V / n,i 

Here the first term is independent of the couplings and is the exact Jacobian due to the 
change of the temporal link variables to L{n) after all the link integrations. The last two 
terms are the leading strong coupling terms in (|) and (^) with the assumption that both 
I and ^ are small and treated on the same footing. Only the leading order terms in f3 
and Pa are retained in the effective action here. Demanding translational invariance for the 
configuration which minimizes the action, one can easily obtain the effective potential for 
the order parameter. Expanding the log term for small L, one has the following results for 
the coefficients, 62, ^4 and fog of the L^, and terms in the effective potential: 

Veff{L) = h2L^ + h^L^ + hL'^ , (15) 
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where 



4 
1 
6 



2 



1 



1 




(17) 



(16) 



(18) 



The positivity of the fee-term (and all other higher terms) ensures that the effective 
potential is bounded from below. Note that for small enough (3 and (3a, ^2 and 64 are also 
positive, favoring thus the confined phase of L = 0. For small Pa and arbitrary /3, 64 remains 
positive but 62 goes through a zero, giving rise to a second order phase transition at a critical 
(3 obtained by setting 62 to zero. As (3a increases, 64 becomes negative above a critical value 
of (3a- The effective potential then has two additional minima in addition to the one at 
L — 0. As (3 increases, these minima deepen and become equal to the one at L = 0, yielding 
a first order deconfinement phase transition. The tricritical point, where the deconfinement 
phase transition changes to become first order, is given by setting the coefficients 62 and 64 
to zero. 

In this leading order strong coupling expansion, the tricritical points are (/3*''»'=^»*^ /3^*'=^**) 
= (0.913,0.649), (1.91,2.418) for Nr=2 and 4 lattices respectively. This suggests that the 
tricritical point moves towards the top right corner {(3 = 00, (3a = 00) of the phase diagram as 
the lattice spacing is reduced. Thus the the tricritical point will not seen by the continuum 
limit. Of course, one needs to improve the leading order strong couphng result as A^t- 
incr eases and check that this conclusion remains unchanged. Nevertheless, these results are 
encouraging for two reasons. Firstly, they provide a concrete example of how the SU{2) 
gauge theory at finite temperature can have a first order phase transition. Indeed, it should 
be noted that the bulk phase transition plays no role above in changing the order of the 
deconfinement phase transition at large (3a- Secondly, the qualitative trend suggested by this 
simple exercise is in agreement with the naive idea of independence of physical results with 
respect to irrelevant couplings. Of course, the key question of the limit of the coefficients 6^, 
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as the lattice spacing a — 0, can only be resolved by simulations at present and, in principle, 
the trajectory of the tricritical point could go either way in that limit. In the next section, 
we describe the results of our simulations which were made in an attempt to answer this 
issue. 



4. RESULTS OF THE SIMULATIONS 



Our Monte Carlo simulations were done using Metropolis algorithm on A^^ x A^,- lattices 
with Nr=2, 4, 6, 8 and Ncr= 8,10,12,16. The many different values of No- were chosen to study 
the finite size scaling behavior of the theory and to compute the critical exponent u, while 
the Nj. values were chosen to monitor the movement of the tricritical point with decreasing 
lattice spacing. We also used N„= 16 in one case to be sure of the critical exponent. The 
possible ranges for the tricritical points for different Nr were known from our earlier work, 
and the simulations were carried out at (3 and i3a in these ranges. Histogramming techniques 
were used to extrapolate to nearby f3 values while estimating the height and location of 
the peak of various susceptibilities. The values of the critical exponents from our earlier 
simulations P^,ITB| on Nt-=2, 4 lattice are summarized in Table 1, where the result of Ref. 



1^ for Pa = 0.0 is also given. One sees a good agreement with the Ising model exponent for 
Nt- = 4 for Pa < 0.9, suggesting that the tricritical point must lie at higher Pa in this case. 
On the other hand, no check of the universality with the three dimensional Ising model has 
so far been made for the Nt-=2 lattices: the corresponding u is unknown although all earlier 
simulations do indicate a continuous transition. From our previous work [|1^ , we know that 



the transition for = 2 and for Pa > 1-4 is a strong first order one. For Pa = 1-4, the 
exponent uj = 3.25(24), with a typical tunneling time of ~ 30,000 — 40,000 Monte Carlo 
sweeps. Thus the effective ranges for the Pa of the tricritical points were 1.4 > for 
Nr = 2, and > 0.9 for Nr = 4. No firm upper bound was known for the latter case, 

although we had good indications that 1.5 > P^J'''^"^, as Ref. [Ill] found a co-existing two 



state signal in both the Polyakov loop, L, and the plaquette, P, at Pa = 1-5. 
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In our earlier work, the simulations at /3a=1-1 on the = 2 lattices did not reveal 
a clear three peak structure in the histogram of the Polyakov loop on the N„ = 8, 10, 12 
lattices, although the peaks did become a bit sharper on going to the N^j = 12 lattice. 
Correspondingly the determination of the critical exponent did not fix the order of the 



transition uniquely. This is similar to the A'^=4 results at /3^=1.1 [11^ where the histograms 
and the evolution graphs of the Polyakov loop gave a very weak two state signal with an 
u ~ 2.34, lying between values characteristic of first and second order phase transition. 
Such a behavior of the deconfinement transition can be understood from the point of view 
of the effective potential in terms of the Polyakov loop, if these simulations were indeed 
close to the tri critical point. As argued in earlier section, the first two leading coefficients 
of and terms of the effective potential are then close to being zero, leading to a 
reasonably flat effective potential around L=0.0, a fact which we will later exploit to conclude 
about the movement of the tricritical point as the lattice spacing is reduced. However, as 
a consequence, much larger statistics is required to sample the exact nature of the effective 
potential near the tricritical point to separate a weak flrst order transition from a second 
order one. We therefore increased the statistics to typically 4 x 10^ sweeps to compute the 
critical exponents and focused more on (3a close to 1.1. We, however, also made simulations 
on Nr=2 lattices at (3a=0-0, 0.8, and 0.9 to determine the critical exponent u and thus the 
range for the tricritical point more precisely. The observables were typically recorded after 
every 20 sweeps to reduce auto-correlation. The errors were estimated by further binning 
the data and the typical bin size was O(IOO). 

Figs. 2a, 2b, 2c and 2d exhibit the results for L-susceptibility for the = 2 lattices for 
(5a = 0.0, 0.8, 1.1 and 1.25 respectively. The results for (3a = 0.9 are similar to those in Fig. 
2c and are therefore not shown here. In each of these flgures, the range of expected values 
for the peak height for the N„=10 and 12 lattices is also shown by two horizontal lines by 
assuming the validity of eq. (|13]), u = 1.97, and by using XNt=8 ^ach case. The errors on 
the respective Xn'^=8 induce the spread between the lines. We always chose a fresh starting 
point in an iterative manner, if the initial guess was too far away from the extrapolated 
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estimate for the location of the peak. This reduced the influence of the unknown systematic 
errors in the peak height due to our extrapolation procedure. As seen in Figs. 2a-2d, we do 
hope that this source of errors has been brought under control by our choice of the simulation 
points, and that it does not annul our conclusions about whether the critical exponent is 
close to 3 or 1.97. 

The independence of uo for (3 a = 0.0 (and small 13a) for =2 and 4 and its agreement 
with 3-d Ising model value is satisfying since there is only one known critical point in Z{2) 
symmetric theories and a priori one expects the exponents to be universal. Moreover, 
strong coupling arguments, which predict a second order phase transition for small /3a, 
should be more reliable for smaller A^^. Quantitatively, however, one notices the leading 
order strong coupling prediction for /5^^**(A^t- = 2)= 0.816 for (3a = 0.0 to be far away from 
the corresponding Monte Carlo determination. Furthermore, its /S^-dependence seems to 
be also in the wrong direction. Thus, one really could have expected surprises in form of 
a qualitative difference from the strong coupling prediction as well. On the other hand, it 
may be more natural to expect the effect of higher orders in quantitative shifts and not in 
qualitative features. Since it is unclear whether iV^=4 is in the strong coupling region, the 
universality conjecture for the critical exponents needs to be tested on lattices with larger 
Nr and thus closer to the continuum limit even for the Wilson action, i.e, (3a = 0.0. 

The values of the critical coupling, (3^^"*' and the finite size scaling exponent, u;, obtained 
by fitting the peak heights to eq. (p!3D , are given in the Table 2 for all the (3a values we 
investigated, including (3a = 0.9. These estimates of uj in Table 2, along with the agreement 
in Fig. 2 with the predictions based on a = 1.97 show the deconfinement phase transition 
for Nt- = 2 lattices to be a clear second order with Ising model exponents for 1.25 > (3a> 0.0 
. At /3a=1-25, additional simulations were performed on A"o-=16 to confirm the second order 
nature of the transition. Note that if the predictions for uj=3, corresponding to a first 
order phase transition were to be displayed in Figs. 2 then they would overshoot by a large 
amount, especially in Fig. 2d where they would be too big by a factor of 2 for the N„ = 16 
lattice. Also interesting to note is the resolution of the ambiguity in establishing the order of 
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the phase transition for (3a = 1.1 in these better statistics simulations. The finite size scahng 
exponent is uj= 1.79(02) and is thus a lot closer to the Ising value. The tricritical point T 
on a Nr=2 lattice is thus constrained to lie in the range 1.4 > (3a > 1.25, as indicated in 
Fig. 1 by the gap between the filled circles for the first order transition points and hollow 
circles for the second order transition points. Thicker dashed and dotted lines show the first 
order and second order deconfinement phase transition lines. 

Since the exponent uo was found [|l4|to be 2.31(28) for = ^ (3a = 1.1 and the 
histograms of the Polyakov loop signaled a very weak first order transition, we chose to 
re-investigate the transition first at /3a = 1.1 and then move to larger values of (3a- Fig. 3 
shows the evolution of both L and P at (3a =1.25 on Na-=8, 10, 12 lattices. These figures 
clearly show the coexistence of two states at the deconfinement transition, since (L) ~ 
for one of the phases while it is nonzero and large for the other. The plaquette, P, has 
a discontinuity at the same location, and further the number of tunnelings and duration 
in each phase do indeed decrease as the spatial volume increases. Figs. 4a and 4b show 
the Polyakov loop susceptibility at /?a=1.1 and 1.25 on A^^^ = 4 lattice for various spatial 
volumes. As in Figs. 2, the expected peak heights for the bigger two lattices are shown by 
horizontal bars. The only difference here is that the solid horizontal lines in Fig. 4-b show 
the predicted Xn'^=io 12 by assuming u = 3.0, while the broken lines in Fig. 4b and the solid 
lines in Fig. 4a are for the Ising value 1.97. The values of couplings at which simulations 
were performed along with the critical values and the fitted u are given in Table 3. The 
most astonishing result is that the transition at Pa =1.25 is a first order transition with 
uj = 3.13(01). This needs to be contrasted with i) the results for N^- = 2 and in particular. 



the Fig. 2d and ii) the results of Ref. p2[ where a clear absence of a first order bulk phase 
transition at Pa =1.25 was shown. 

These simulations thus indicate that the tricritical point for A'^^= 4 lies definitely below 
(3a = 1.25 whereas the corresponding Nr = 2 point is definitely above (3a = 1.25. This is 
also clearly seen in Fig. 1, where again the first order transition points for A',- = 4 are shown 
by filled squares, the second order points by hollow squares and the gap between them is 
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the allowed range for the tricritical point. As one can see, the tricritical point does shift as 
the temporal lattice size increases from 2 to 4. However, the direction of the shift is almost 
orthogonal to the strong coupling prediction of the previous section and its magnitude is 
also much smaller. The shift, on the other hand, suggests a possible lack of any correlation 
of the bulk transitions, if any, with the key observation of the change of the order of the 
deconfinement phase transition. This is so since any possible bulk transition for (3a < 1.25 is 
definitely not [^] a first order phase transition, making it an unlikely cause of the behavior 
seen in Figs. 3 and 4. A plausible explanation then is indeed the phenomenon seen in 
the strong coupling calculations in Sec. 3, i.e., a change of sign in the coefficient 64 of the 
effective potential for the order parameter, L. 

The above downward movement of T in the extended coupling plane is puzzling. It calls 
for a more detailed cross-check on the 3-d Ising model universality, especially for positive 
values of 13a but also perhaps for larger A',- at (3a = 0.0. Note that the latter has so far 
been demonstrated to a very good accuracy only on the A^,- = 2,4 lattices. Therefore, we 
decided to monitor the deconfinement transition further by simulating the model on A^,- = 6 
and Nt-=8 lattices with A^o-=12 and 16 respectively. On these lattices the simulations were 
performed at /?a=1-1 and 1.25 to determine the range in which T may lie. At /9^=1.25 on 
A^^ =6 lattice we found the transition to be first order. The corresponding histogram is 
plotted in Fig. 5 at (3=1.2184. The three peaks are clearly visible and distinct, though not 
of equal height. This figure suggests the transition point to be at a slightly higher value 
of (3 than 1.2184. In choosing this (3 we were guided by the location of the peak of the 
L-susceptibility to locate the critical point. We have found that this criterion to determine 
the critical point differed a little from that of the effective potential picture (i.e, the nature 
of the L-histograms) , although they will coincide in the thermodynamic limit. Comparing 
the shapes of the histograms at (3a=1-1 in Figs. 6-a,b,c (A^t-=4,6,8 respectively) on the other 
hand, one observes that their profiles at the peak tend towards a Gaussian behavior as one 
makes the lattice spacing smaller by going from A^ = 4 to 6 to 8. This is the expected 
behavior if the tricritical point shifts upwards from (3a — 1-1 for A^,- = 4. Since the effective 
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potential (|1^) will have 62 ~ &4 ~ at the tri critical point, it will have a reasonably flat 
bottom, causing the very flat top for the Nt-=A histogram. In Fig. 6-a we have plotted the 
histograms for A^^=4 and Ncr=8,10, 12 to conflrm the closeness of the tricritical point with 
the run point, i.e. /3=1.327 and 1.32685. Moreover, one clearly sees that the flatness is not 
a flnite N„ artifact. The reduction in the above flatness would signal 64 becoming nonzero, 
as in a typical second order transition. More quantitatively, the fluctuations in the Polyakov 
loop at (3c = 1-327 on the A^,- =4 lattice, are in the range AL ^ ± 0.3, while AL ^ ± 0.1 
on the Nt = 6 lattice at /5c=1.339, with a reduction by a factor of 3 in the flatness of the 
histogram at the top. The statistics for the above runs on A^,- =4 and 6 lattices was roughly 
2 X 10^ and 7 x 10^ respectively. In principle, this could be taken as a hint that the tricritical 
point has moved up. However, one knows that L measures the free energy of a point-like 
test charge and thus has A^,--dependent corrections. These reduce the value of L just above 
the phase transition as increases. A perturbative estimate of these corrections will 
presumably still not yield the factor of two but unless it is quantitatively so demonstrated 
one is handicapped in drawing a flrm conclusion about lack of the flnite A^t- effects in the 
flatness of the effective potential. The Polyakov loop histogram in Fig. 6-c corresponding 
to A^T- = 8 has no indication of the co-existence of 2 phases and therefore no evidence of a 
tricritical point in its vicinity. In this flgure we have plotted 2 histograms corresponding to 
(3=1.35 and 1.3508 showing the the sharp behavior of the transition. This indicates that the 
qualitative behavior of the histogram changes drastically with the coupling. This again is not 
expected near the tricritical region. The above drastic reduction in the flat region compared 
to Nt =4 lattice could again be partially due to the same cause as above. Nevertheless, 
these qualitative observations suggest that that the tricritical point T moves upwards as Nr 
is increased from 4 to 6 and then from 6 to 8. Even assuming that this upward motion of 
the the tricritical point continues, one will need a lot larger lattices to conflrm universality 
for (3a > since the upward movement is rather small. 

5. SUMMARY AND DISCUSSION 
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We simulated the extended action of eq.([l|) on A^,- = 2, 4, 6 and 8 lattices with varying 
spatial sizes and determined the order of the deconfinement phase transition by obtaining 
the critical exponent uj using finite size scaling theory and also by observing the behavior 
of the histograms. For (3a = 0.0 and N-r = 2, we find uj = 1.92 ± 0.01 which is in good 
agreement with the corresponding Ising model value 1.97 ± 0.03. Thus the deconfinement 
phase transition for the Wilson action for both = 2 and 4 is of second order and is in 
the same universality class as the three-dimensional Ising model. Our simulations show this 
universality to persist when (3a is turned on. For the = 2 lattice, this is true for at least 
(3a < 1-25 while for Nr = 4 the phase transition at (3a = 1-25 is already of first order, with 
C(j = 3.13 ± 0.01. The tricritical point where the deconfinement phase transition changes it 
order is therefore definitely above (3a = 1-25 for = 2 while it is definitely below (3a = 1-25 
for Nt- = A lattices. Placing it in the middle of the ranges we have narrowed down, it shifts 
vertically downwards by about 0.15. There is a horizontal positive shift in (3, of the order 
of 0.09, as well. 

The comparison of the above shift of the tricritical point on Nt-=2 and 4 lattices with 
the leading order strong coupling prediction shows surprises. It is, of course, reassuring that 
the change of the order and the existence of tricritical points for each N^- is as predicted. 
However, the predictions do very poorly on a quantitative level. In particular, the directions 
of the shift of the phase transition as (3a is turned on, and more importantly, the predicted 
vertically upward shift of the critical point are in complete contrast with the simulation 
results for A^^- =2 and 4. On the other hand, the leading order strong coupling equations 
are known to fail quantitatively for the Wilson theory for both A^^- = 2 and 4 as well. 
Putting (3a=0-0 in eq. (|1^), we recover p8[ the criticality condition for the Wilson action 



: = 4(1/24)^. Therefore, the values of the critical couplings to this order on 

Nt- = 2 and 4 lattices are 0.816 and 1.807 respectively. These should be compared with 
the corresponding Monte Carlo values of this work {(3c{Nr = 2) = 1.88(01)) and of Ref. 
|T^(/9c(A^T = 4) = 2.30(01)). It would be interesting to check whether the inclusion of higher 



orders in the strong coupling expansion can yield a better agreement with the simulation 
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results, especially for the direction of the vertical shift of the tri critical point. 

The above downward movement of the tricritical point on going from Nt-=2 to 4 was 
also observed in Ref. where the Villain form of eq. (|T]) was simulated. This unexpected 
behavior of the tricritical point may therefore need to be taken seriously. If this trend 
persists even on larger lattices then the continuum limit of SU(2) lattice gauge theory could 
correspond to a first order deconfinement transition. Hence we simulated this model on larger 
lattices with A^,- =6 and 8 and found some hints of an upward movement of the tricritical 
point by comparing the shapes of the Polyakov loop histograms. While it unfortunately 
is not very conclusive, it is encouraging that the trajectory of the tricritical point could 
potentially be turning up in the coupling plane for these A^^. The tricritical point for 
both these lattices was still found to be below (3^=1.25 and above (3a = l-l- Much more 
computational resources on bigger lattices with A^,- much bigger than 8 are necessary to 
confirm this. 

It may therefore be important to understand and explain the origin of the change in 
the order of the deconfinement transition even away from the continuum, especially since 
the results of Ref. p2] suggest a lack of a first order bulk phase transition at (3a = 1-25 



for this action and thus make it implausible that a bulk transition is responsible for such a 
change. To answer this question, we plan to consider the Villain form for the SO (3) part of 
the extended action @]: 

Z= E [UdU,{n)exp[J2((^ + (^a,]TrpU, + Xa,] , (19) 

Here dp is a Z2 plaquette field and the summation over it ensures the invariance of the second 
term above under U^{n) —U^{n). For A = 0, this action is again in the same universality 
class as that of (|I|). In fact as already mentioned, its simulations ||T6[ on A^,- = 2,4 lattices 
led to exactly the same behavior of the tricritical point as reported in this paper. Besides 
computational advantages [|r^, the theoretical advantage of this action is that unlike eq.(0), 
the SO (3) monopoles and their dynamics is manifest in the form of the Z2 plaquette field. 
The SO (3) monopole charge density is given by HpGc c"pj here the product is over the 6 faces of 
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a cube 0. The last term in the above equation is the potential energy for these topological 
degrees of freedom. In the extreme (A oo) case when all the 50(3) monopoles are 
suppressed (VcTp = +1), the above extended action reduces to Wilson action with redefined 
coupling and therefore has only a second order deconfining transition on small A^^ lattices. 
Therefore, these topological degrees of freedom may have a crucial role in changing the order 
of the transition. In the extended coupling plane these monopoles cost less and less energy 
as the adjoint coupling is increased with decreasing values of /3. Therefore, above 
they might condense giving rise to a first order transition. This can be checked by simulating 
the above model. If true, this phenomenon will be particularly appealing because precisely 
the same degrees of freedom and their condensation have been attributed to the first order 
nature of the 'bulk transition' which we find to be first order deconfinement transition. It 
may thus also enable in resolving the physical nature of the transition. 
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TABLES 

TABLE L The average values of the critical exponent lo at different values of the adjoint 
couplings. The expected value is 1.97 (3.0) if the deconfining phase transition is second order (first 



order). The data taken from Refs. |12,14|,1 



13a 


Nr 


UJ 


0.0 


4 


1.93(03) 


0.5 


4 


1.92(29) 


0.75 


4 


1.53(32 


0.9 


4 


2.10(22) 


1.1 


4 


2.34(15) 


1.4 


2 


3.246(243) 



TABLE II. The values of (/?, Pa) on x 2 lattice at which simulations were performed, P'^"*- 
and the finite size scaling exponent lo. The expected value for uj is 1.97 (3.0) if the deconfining 
phase transition is second order (first order). 



Pa 



P 



Pc 



0.0 



10 
12 



1.90 
1.885 
1.877 



1.88 
1.878 
1.877 



1.92(01) 



0.8 



10 
12 



1.368 
1.367 
1.368 



1.368 

1.3664 

1.366 



2.03(01) 



0.9 



10 
12 



1.3 
1.3 
1.3092 



1.31 
1.31 
1.3088 



1.83(02) 



1.1 



10 
12 



1.201 
1.2 
1.1999 



1.201 
1.2 
1.1995 



1.79(02) 



1.25 



10 
12 
16 



1.12 
1.12 
1.12 
1.12 



1.1203 
1.1203 

1.12 
1.1196 



1.97(01) 



TABLE III. The values of {(3, Pa) on x 4 lattice at which simulations were performed, 
and the finite size scaling exponent u. 



8 1.327 1.327 

1.1 10 1.327 1.3274 2.11(02) 

12 1.32715 1.32715 

8 1.2146 1.2142 

1.25 10 1.214 1.2144 3.13(01) 

12 1.2144 1.2143 



FIGURE CAPTIONS 



Fig. 1 The phase diagram of the extended SU(2) lattice gauge theory. The sohd hnes are 
from simulations done on a 5^ lattice by Bhanot and Creutz [Q]. The light dashed line 
indicates the absence of first order bulk transition below (3 a = 1.25. The dotted(thick 
dashed) lines with hollow(filled) symbols are the second(first) order deconfinement phase 
transition lines on A^^ =2 (circles) and 4 (squares) lattices. 

Fig. 2 Polyakov loop susceptibility at (a) (3a=0.0, (b) (3a=0.8, (c) ^a=1-1 and (d) /5a=1.25 
on 8^ X 2, 10^ X 2 and 12^ x 2 lattices. At I3a = 1.25 result on 16^ x 2 lattice is also shown. 
The points with error bars are results of simulations and the curves are extrapolations by 
the histogramming technique. The horizontal lines are predictions assuming a second order 
deconfinement transition, as explained in the text. 

Fig. 3 Evolution of |L| and P at /3a = 1.25 on (a) 12^x4 (/5 = 1.2144), (b)103x4 (/3 = 1.214) 
and (c)83 x A {/3 = 1.2146) lattices. 

Fig. 4 Polyakov loop susceptibility at (a)/?A=l.l and (b)/3A=1.25 on 8^ x 4, 10^ x 4 and 
12^ X 4 lattices. The solid horizontal lines are predictions assuming a first (second) order 
deconfinement transition at Pa = 1.25 (/9a=1.1), as explained in the text. The broken lines 
in Fig. 4a are predictions for second order. 

Fig. 5 

The probability density of L at Pa = 1-25 on 12^ x 6 lattice at /3=1.2184. 
Fig. 6 

The probability density of L at (3a = 1.1 on x Nr lattice for a] A^=4, A^=8,10,12 
and /5=1.327,1.327 and 1.32685 respectively, b] A^=6, = 12 and /3= 1.339, c] A^=8, 
= 16 and p= 1.35 and 1.3508. 
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